Extraction of the homogeneous linewidth of a fast spectrally diffusing line 
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We present a simple method to extract the homogeneous linewidth of a single photon emitter line 
exhibiting fast (down to 1 ns) spectral diffusion (SD). It is based on a recently developed technique 
using photon correlation measurements on half of the line. Here we show that the SD induced 
bunching depends on the ratio between the width of the homogeneous line and the spectral diffusion 
amplitude. Using this technique on a CdSe/ZnSe quantum dot, we investigate the temperature 
dependence of its fast SD amplitude and its homogeneous excitonic linewidth. 

PACS numbers: 78.67.Hc, 78.67.Uh, 78.55. Et, 42.50.Lc, 



Spectral diffusion (SD) corresponds to random spec- 
tral jumps of a narrow line within a broader spectral 
profile. It is a prominent issue in spectroscopy [T}{5]. 
Prom a practical point of view, it broadens the observed 
linewidth and prevents the access to the intrinsic line 
properties. On a more fundamental side, it gives some 
very local information on the microscopic environment of 
a single light emitter embedded in a solid matrix, or mov- 
ing within a fluid or a gas. The quest for large light mat- 
ter coupling in condensed matter physics makes use of 
photonic nanostructures such as photonic crystals [5lll0|. 
photonic wires |lllH^] or metallic nanoantenae [HI] where 
the emitter is not far from a surface and may therefore 
undergo spectral diffusion caused by surface charge fluc- 
tuations. 

Visualizing directly the spectral wandering by record- 
ing a time series of spectra has been so far the usual 
method to observe SD [2-8 . For single photon emitters, 
the time resolution was therefore limited by the mini- 
mum time of about 1 ms required for a photon counting 
charged coupled device (CCD) to acquire a spectrum. 

Except from the exotic possibility of statistics such as 
Levy flight [T3], SD is usually entirely described by its 
fluctuation amplitude and its characteristic time. But, 
additionally, one might wish to access the width of the 
homogeneous wandering line. When the characteristic 
SD time is faster than the time necessary to record a spec- 
trum, this information is hidden by the inhomogeneous 
linewidth and difficult to access for a single photon emit- 
ter. Sophisticated methods featuring both high spectral 
and temporal resolution are required to access this quan- 
tity. To our knowledge, photon correlation Fourier spec- 
troscopy (PCFS) [15 is the only published method able 
to extract the homogeneous linewidth of a rapidly wan- 
dering line of a single emitter together with its SD param- 
eters. In PCFS, high spectral resolution (AA/A = 10~ 6 ) 
is provided by a Michelson interferometer, and temporal 
resolution of 20 /is is achieved with a photon correlation 
set-up at the outputs of the interferometer. Alternatively, 
four wave mixing |16| is able to give the coherence time of 
the transition of a single emitter undergoing spectral dif- 



fusion, but without any access to SD characteristic time. 
In the case of an ensemble of emitters, lock-in hole burn- 
ing techniques [17J can give the homogeneous linewidth. 

Our recently developed photon correlation spec- 
troscopy (PCS) technique [18] converts spectral fluctu- 
ations into intensity fluctuations, as also reported in 
[151 Hi?l |2"U] . This simple method benefits from the sub- 
nanosecond time resolution of an Hanbury-Brown and 
Twiss photon correlation set-up and the spectral resolu- 
tion of the spectrometer. In short, it is based on correla- 
tions of photons emitted within a spectral window nar- 
rower than the SD broadened line. Owing to the wander- 
ing of the homogeneous line, the emission energy stays a 
limited time within this spectral window leading to pho- 
ton bunching (see fig [T| . The characteristic time r c of 
this effect can be easily accessed by photon correlation 
even though it is much shorter than the inverse photon 
count rate. In practice, spectral diffusion in the 1 ns 
range can be measured even though the count rate is less 
than one photon every 10 fj,s. This opens a new domain 
of fast spectral diffusion phenomena that was previously 
not accessible. 

In this letter, we present a development of the PCS 
technique that enables the extraction of the homogeneous 
linewidth of the spectrally diffusing line of a single pho- 
ton emitter, together with its complete SD parameters. 
The method is performed on CdSe quantum dots (QDs) 
embedded in a ZnSe nanowire (NW) [21 , with a single 
object photoluminescence linewidth between 1 and 2 meV 
|22| . After describing the principle of our measurement, 
we discuss the properties of the half-line autocorrelation 
function (HLAF) of a spectrally diffusing single photon 
emitter with a finite homogeneous linewidth. We show 
that the contrast of the measured photon bunching is 
related to the ratio between the SD amplitude and the 
homogeneous linewidth of the wandering line. We evalu- 
ate these two quantities separately on our NW-QD and 
study their behavior when temperature is increased as 
an illustration of the possibilities offered by our method. 
This allows us to probe the emitter energy statistic as a 
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Figure 1: (Color online) Principle of the PCS method in the 
case of infinitely narrow line (a-c) and in the case of a line 
of finite width (d-e). In a) and d) are represented the spec- 
tra of a spectrally diffusing single emitter with an infinitely 
narrow linewidth (in a)) and a finite linewidth (in d)). Only 
the right (R) spectral window is selected for detection. This 
corresponds to the photon statistics of the R side as shown 
in b) and e) and to the half line autocorrelation functions 
displayed as a solid line in c) and f). In b) and e), "red pho- 
tons" correspond to photons emitted from a line whose center 
is located on the R side. However, in e), owing to the finite 
linewidth of the homogeneous line, red (resp. blue) photons 
can be emitted on the L (resp. R) side, leading to a reduction 
of the SD induced bunching. In c) and f), the zero delay dip 
is due to the single photon type antibunching and the dotted 
line is the bunching contribution caused by spectral diffusion 
only. 



function of temperature. 

In the case of an infinitely narrow homogeneous line 
as shown in fie (l|a-c), the HLAF, noted g^(R,r) when 
the R side of the inhomogeneous line is detected, writes 



7 (2) 



(R,t) = [l + /3exp(- 7c r)][l-exp(-(r + 7)r)], (1) 



with 7 C = — being the SD correlation rate, 7 the ra- 
diative emission rate, r the pumping rate and /? the SD 
induced bunching factor. This expression is the product 
of two terms. The last term corresponds to the anti- 
bunching of the single photon emitter. The first term 
accounts for spectral diffusion and its influence on the 
ability of the emitter to send photons in the R spectral 
window. Note that /3 = 1 in the case of an infinitely 
narrow homogeneous line when the two energy sides are 
of equal importance (ie half of the total count in the R 



side). Equal distribution between L and R is checked 
very carefully since deviation from this situation leads to 
a different value of /3 |18j . 

We now consider the case of a finite homogeneous 
linewidth (fig(ljd-f)). In this situation, even if the ho- 
mogeneous line is centered in the L side, a photon can 
nevertheless be detected in the R side and vice versa as 
represented in fig. [lje). This randomizes the photon 
stream and results in a reduction of the bunching con- 
trast f3 as seen in figjljf). 

A probabilistic model was developed considering 
the consequence of a finite linewidth on the HLAF. The 
key idea is to take into account the combined effects of 
the correlated energy statistics of the SD together with 
the uncorrelated energy statistics introduced by the ho- 
mogeneous linewidth. When the photons are detected on 
the R side, the bunching contrast f3 is then found as 



{oti 



OLout ) 



(2) 



\pt>in ^out ) 

where ai nmOU t is a normalized overlap integral defined as 
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where D^{ii) is the SD line distribution of full width 
at half maximum (FWHM) S, i.e. the probability den- 
sity associated with the position of the homogeneous line 
median energy fi, and H a (E — li) the homogeneous line 
shape with median li and FWHM a. We assume here 
that each of these function has a given shape that can 
be described by a single width parameter (£ for Dy, and 
a for H a ). The coefficients di n and a ou t are the proba- 
bility to detect a photon in the R spectral window with 
the homogeneous linewidth respectively centered in R or 
in L. They only depend on the homogeneous and SD line 
profiles H a and Dy,. If the a priori knowledge of the 
physical system allows us to chose the shape of the SD 
distribution and of the homogeneous line, which is gen- 
erally the case, the bunching contrast f3 will only depend 
on their linewidth ratio er/S. This is the main result of 
the paper. 

In the following, we have used a Lorentzian shape of 
FWHM a for H a and a Gaussian shape of FWHM £ for 
Dy, to model the bunching factor j3 and fit the experimen- 
tal data. We are aware that a Lorentzian shape is only an 
approximation of the real homogeneous line shape, espe- 
cially at elevated temperature when the exciton-phonon 
coupling leads to more complex line shapes but this 
gives nevertheless a meaningful estimation of the homo- 
geneous linewidth and its variation. The Gaussian shape 
of the spectral diffusion distribution is justified by the 
Kubo- Anderson theory which describes the spectral dif- 
fusion as the result of fluctuations of a very large number 
of independent identical random variables, leading glob- 
ally to a Gaussian distribution of the emitter energies in 
the "slow fluctuations" regime where Er c 3> ft [25-27]. In 
our case £ sa 1 meV and t c ss 10 ns so that £r c sa 10 4 ft. 
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Figure 2: (Color online) The blue solid line is the theoretical 
value of the bunching factor /3 (cf eq |2|) as a function of the 
linewidths ratio cr/E in the case of a Gaussian inhomogeneous 
line with FWHM E and a Lorentzian homogeneous line with 
FWHM a. This result is confirmed by a Monte-Carlo simu- 
lation (red dots). The green dashed line is the relative error 
<5((j/E)/((j/E) for the actual experimental error 6(3 = 0.04. 



Within the above assumptions, the bunching factor /3, 
computed from eq.([2]), is plotted in fig (2] as a function 
of cr/E together with a Monte-Carlo simulation. The 
bunching factor decreases with this ratio because the 
uncorrelated statistics introduced by the homogeneous 
linewidth is degrading the bunching induced by the cor- 
related SD process. Note that for an infinitely small ho- 
mogeneous linewidth the bunching factor is equal to unity 
(a in = 1/2 and a out — 0), the energy of the emitter being 
only defined by SD . On the contrary, for an homogeneous 
linewidth larger than the SD amplitude the bunching 
vanishes (cvi„ — > 1/4 and a ou t — > 1/4), as the uncorre- 
lated homogeneous linewidth process dominates and the 
emitter energy statistics in the time domain tends to be 
Poissonian. As seen with the green dashed line in figj2] 
a relative error on cr/E better than 30% is obtained for 
0.04 < cr/E < 1.1 in our experimental conditions. 

The measurement of the HLAF gives access to the 
bunching factor and therefore to the cr/E ratio. The 
spectrum also depends on these two parameters as it is 
a convolution of the homogeneous linewidth and of the 
spectral diffusion distribution. It leads to a Voigt pro- 
file when the homogeneous and the inhomogeneous line 
are respectively Lorentzian and Gaussian as mentionned 
above. Thus, the recording of the spectrum and of the 
HLAF (cf fig{3| allows the extraction of the homogeneous 
linewidth a together with the SD amplitude E and char- 
acteristic time t c in the 10 ns range. 

Our experimental results have been obtained with a 
time resolved microphotoluminescence setup. The ex- 
citation source is a continuous wave laser emitting at 
A = 405 nm (ie 3.06 eV). The QD luminescence is de- 
tected through a SX — 0.06 nm (ie 5E = 300 fJ,eV) reso- 
lution spectrometer on a charged coupled device (CCD) 
camera for spectra or on a Handbury-Brown and Twiss 
photon correlation set-up of overall time resolution of 800 
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Figure 3: (Color online) (a) Typical measured HLAF at 4K 
together with its best fit. The limited depth of the central 
antibunching dip is due to the limited time resolution of the 
photon correlation set-up. (b) Emission line spectrum fitted 
with a Voigt profile. 



ps (FWHM). 

The measurement method was applied on three differ- 
ent QDs with different SD amplitudes and the extracted 
homogeneous linewidths were around 300 [ieV for the 
three measurements. This actually corresponds to the 
experimental spectral resolution imposed by the spec- 
trometer. This is the main limitation of this method. 
However, we underline the fact that the measured cor- 
relation times are of the order of 10 ns and that it is 
the only reported method able to probe experimentally 
the homogeneous linewidth of such a fast diffusing single 
photon emitter. 

When the temperature is increased, the first noticeable 
consequence on the HLAF is the decrease of the bunching 
contrast (figjlja)). The first interpretation of this obser- 
vation is the increase of the homogeneous linewidth ow- 
ing to temperature induced phonon broadening [231 12H] ■ 
This process becomes dominant above 60 K for this QD, 
and the bunching collapses completely for higher tem- 
peratures meaning that the cr/E ratio becomes large. At 
high temperature the energy of the emitter is less and less 
time correlated and tends to adopt a Poissonian statis- 
tics. We display the extracted homogeneous linewidths 
and spectral diffusion amplitudes (figElb)) together with 
the spectral diffusion rate 7 C — 1/t c (ngj4]c)) as a func- 
tion of temperature for two different QDs with different E 
(blue: QD1, E = 1.7 mcV and red: QD2, E = 1.2meV). 
The effective homogeneous linewidths for both QDs in- 
crease with temperature. They reach the same order of 
magnitude as the fluctuation amplitude for temperatures 
around 40-50 K. For temperatures higher than 60 K, the 
bunching becomes no longer measurable and the cr/E ra- 
tio can no longer be accurately extracted. 
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A remarkable feature that appears in fig|4]b) is that the 
SD amplitude is constant with temperature for the two 
investigated quantum dots. So our results show that the 
main contribution to the QD linewidth is SD for temper- 
atures below 60 K and phonon broadening above. This 
relatively surprising behaviour can be explained with 
the help of the Kubo Anderson model adapted to non 
symmetrical fluctuations |27l 125] . Even though the de- 
scription of spectral diffusion mechanisms in QDs is not 
the main focus of this paper, we give here a brief in- 
terpretation of this result. Emission energy fluctuations 
in semiconductors are caused by captures and escapes 
of carriers in defects in the QD vicinity. These car- 
rier motions are mainly governed by phonon and Auger 
effects. The Auger contribution is similar for capture 
and escape processes and is dominant at high pump- 
ing power. This means that, under our experimental 
conditions, 7^ sa 7^ whatever the temperature. Since 
the SD amplitude is governed by 7^/7^ [27 , it should 
therefore be temperature independent, as observed ex- 
perimentally. Furthermore this effect is rather robust to 
a departure of the 7^/7^ ratio from unity owing to the 
fact that the SD amplitude features a rather flat max- 
imum around 7-1-/74. = 1. Figj4] c) shows that the SD 
rate 7 C = l/r c increases with temperature. This rate is 
given by 7 C = 7^ + Tf [27 . Its rise with temperature 



is caused by the rise of 7^ and 7-f with the increasing 
phonon population. The difference between the two QDs 
can be attributed to a difference in the Auger rate and/or 
to a different defect energy depth. 

As a conclusion, we have demonstrated that the mea- 
surement of the HLAF brings informations on the energy 
statistics of the emitter in addition of providing a sub- 
nanosecond resolution access to the correlation time of 
the fluctuations. The main result of this paper is the 
possibility of extracting of the homogeneous linewidth of 
a fast diffusing line with a simple photon correlation set- 
up. This work brings a new tool for the understanding 
of single emitters environment fluctuations. As an illus- 
tration, we showed that, in CdSe/ZnSe NW QDs, SD 
amplitude does not depend on temperature and that the 
line broadening is only due to the coupling of the emitter 
with phonons. This contributes to a better interpretation 
of the temperature behaviour of QD linewidths. 

We acknowledge support from the French National 
Research Agency (ANR) through the Nanoscience and 
Nanotechnology Program (Project BONAFO ANR-08- 
NANO-031-01) that provided a research fellowship for 
MdH. ME J acknowledges financial support from the 
Nanosciences Foundation "Nanosciences, aux limites de 
la nanoelectronique" (RTRA). 



J. R. Klauder and P. W. Anderson, Phys. Rev. 125, 912 
(1962) 

W. P. Ambrose, and W. E. Moerner, Nature 349, 225 
(1991). 

S. A. Empedocles, D. J. Norris, and M. G. Bawendi, 
Phys. Rev. Lett. 77, 3873 (1996). 

H. D. Robinson and B. B. Goldberg, Phys. Rev. B 61, 
R5086 (2000). 

T. Plakhotnik and D. Walser, Phys. Rev. Lett. 80, 4064 
(1998) 

J. Seufert, R. Weigand, G. Bacher, T. Kummell, A. 
Forchel, K. Leonardi, and D. Hommel Appl. Phys. Lett. 
76, 1872 (2000) 

V. Tiirck, S. Rodt, O. Stier, R. Heitz, R. Engelhardt, U. 
W. Pohl, D. Bimberg, and R. Steingruber, Phys. Rev. B 
61, 9944 (2000). 

L. Besombes, K. Kheng, L. Marsal, and H. Mariette, 
Phys. Rev. B 65, 121314 (2002). 

D. Englund, A. Faraon, I. Fushman, N. Stoltz, P. Petroff 
and J. Vuckovic, Nature (London) 450, 857 (2007) 
A. Majumdar, E. D. Kim, and J. Vuckovic, Phys. Rev. 
B 84, 195304 (2011) 

J. Claudon, J. Bleuse, N.S. Malik, M. Bazin, P. Jaf- 
frennou, N. Gregersen, C. Sauvan, P. Lalanne and J.-M. 
Gerard, Nature Photon. 4, 174 (2010). 

I. Yeo, N.S. Malik, M. Munsch, E. Dupuy, J. Bleuse, Y.- 
M. Niquet, J.-M. Gerard, J. Claudon, E. Wagner, S. Sei- 
delin, A. Auffeves, J. -Ph. Poizat, and G. Nogues, Appl. 
Phys. Lett. 99, 233106 (2011) 

L. Novotny and N. van Hulst, Nature Photon. 5, 83 
(2011) 



[14] X. Brokmann, J. -P. Hermier, G. Messin, P. Desbiolles, J.- 
P. Bouchaud, and M. Dahan, Phys. Rev. Lett. 90, 120601 
(2003) 

[15] X. Brokmann, M. Bawendi, L. Coolen, J. -P. Hermier, 
Opt. Express, 14, 6333 (2006); L. Coolen, X. Brokmann, 
and J.-P. Hermier, Phys. Rev. A 76, 033824 (2007); L. 
Coolen, X. Brokmann, P. Spinicelli, and J.-P. Hermier 
Phys. Rev. Lett. 100, 027403 (2008). 

[16] B. Patton, W. Langbein, U. Woggon, L. Maingault and 
H. Mariette, Phys. Rev. B 73, 235354 (2006); W. Lang- 
bein and B. Patton, J. Phys. Condens. Matter 19 295203 
(2007); W. Langbein, Rivista Nuovo Cimento 33, 255 
(2010) 

[17] P. Palinginis, S. Tavenner, M. Lonergan, and H. Wang 

Phys. Rev. B 67, 201307 (2003). 
[18] G. Sallen, A. Tribu, T. Aichele, R. Andre, L. Besombes, 

C. Bougerol, M. Richard, S. Tatarenko, K. Kheng, and 

J. Ph. Poizat, Nature Photon. 4, 696 (2010) 
[19] A. Zumbusch, L. Fleury, R. Brown, J. Bernard, and M. 

Orrit, Phys. Rev. Lett. 70, 3584 (1993). 
[20] L. F. Marshall, J. Cui, X. Brokmann, and M. G. Bawendi, 

Phys. Rev. Lett. 105, 053005 (2010). 
[21] M. Den Hertog, M. Elouneg-Jamroz, E. Bellet-Amalric, 

S. Bounouar, C. Bougerol, R. Andre, Y. Genuist, J. -Ph. 

Poizat, K. Kheng and S. Tatarenko, J. Appl. Phys. 110, 

034318 (2011). 

[22] S. Bounouar, C. Morchutt, M. Elouneg-Jamroz, L. Be- 
sombes, R. Andre, E. Bellet-Amalric, C. Bougerol, M. 
Den Hertog, K. Kheng, S. Tatarenko, and J. Ph. Poizat, 
Phys. Rev. B 85, 035428 (2012). 



[23] A. Trichet and S. Bounouar. arXiv: 1201.4016 



5 



■),. 



CC 5 -20 3 10 



— ~~"A / — ~~ 



Delay r (ns) 



b) 



£1 

B 
W 



-ile-***-*^^ 

aoo 



c) 
v"0.2 

CO 
fZ 

"io.o. 



O §* 



+ 4 , + : + 

20 40 60 
Temperature (K) 



[24] L. Besombes, K. Kheng, L. Marsal, and H. Mariette 

Phys. Rev. B 63, 155307 (2001). 
[25] R. Kubo, J. Phys. Soc. Jap. 9, 935 (1954). 
[26] P. W. Anderson, J. Phys. Soc. Jpn. 9, 316 (1954). 
[27] A. Berthelot, I. Favero, G. Cassabois, C. Voisin, C. De- 

lalande, Ph. Roussignol, R. Ferreira, and J.-M. Gerard, 

Nature Physics 2, 759 (2006). 
[28] I. Favero, A. Berthelot, G. Cassabois, C. Voisin, C. De- 

lalande, Ph. Roussignol, R. Ferreira, and J.-M. Gerard, 

Phys. Rev. B 75, 073308 (2007). 
[29] R. J. Cook and H. J. Kimble, Phys. Rev. Lett, 54, 1023 

(1985). 



Figure 4: (Color online) a) HLAF for different temperatures 
for QD1 with the corresponding fit. b) Temperature depen- 
dency of the homogeneous linewidth (o) and SD amplitude (*) 
extracted from HLAF and linewidth measurements for QD1 
(in blue) and QD2 (in red), c) Temperature dependency of 
the SD rate >y c for QD1 (in blue) and QD2 (in red). 



